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Abstract 

At the cohiders RHIC and LHC, nuclei at the ultrarelativistic energies of 100 GeV/A 
and 2.7 TeV/A will be smashed together with the hope of creating an elusive and 
short-lived state of matter called the quark gluon plasma. The initial conditions 
which determine the dynamical evolution of the quark gluon matter formed in the 
central region after the collision depend crucially on the small x component of the 
nuclear wavefunction before the collision. In this comment, we discuss recent work 
which argues that, for large nuclei, weak coupling techniques in QCD can be used 
to calculate the distribution of these small x, or wee, partons. The ramifications of 
this approach for the dynamics of heavy ion collisions and the various signatures of 
a quark gluon phase of matter are discussed. 



1 Introduction 



What does a nucleus look like when it is boosted to ultrarelativistic energies? The 
special theory of relativity tells us that the nucleus must contract a distance -R/7 in 
direction of its motion, where R is its radius and 7 >> 1 is the Lorentz factor. If we 
increase 7 indefinitely, do we expect the longitudinal size of the nucleus to shrink to 
a point? What does this statement mean in terms of the underlying parton degrees 
of freedom? What happpens to its transverse size-does it approach a constant at 
asymptotic energies or does it keep growing [j^? 

With the advent of the Relativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider (LHC), which will collide large nuclei at ultrarelativistic energies, 
the above questions are not merely academic but are extremely relevant to under- 
standing the initial conditions for these collisions. The primary objective of heavy 
ion collisions at these energies is to investigate the possible formation of a soup of 
quark gluon matter-often simply called the quark gluon plasma-and a phase transi- 
tion of the plasma to hadronic matter . The formation of the plasma and indeed 
the dynamics of any subsequent phase transition to hadronic matter, will depend 
sensitively on these initial conditions. 

In this comment, I will discuss recent work p|-[lO[| which seeks to answer the 
above questions quantitatively by addressing the problem of initial conditions for 
nuclear collisions within the the framework of Quantum Chromodynamics (QCD). 
The center of mass energies of the colliding nuclei at RHIC and LHC are 100 GeV 
and 2.7 TeV respectively. Since these energies are far greater than the typical 
energies of nuclear interactions, the appropriate degrees of freedom in describing 
these relativistic nuclei must be quarks and gluons, whose interactions are described 
by QCD. 

For ultrarelativistic nuclear collisions at central rapidities, the properties of 
quarks and gluons at very low values of a; ~ kt/^/s are relevant. (Note that x is the 
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light cone momentum fraction of the nuclear momentum carried by the quark or 
gluon, kt is its transverse momentum and ^/s is the center of mass energy). Recently, 
there has been renewed interest in QCD at small x because of the results of the 
deeply inelastic electron proton scattering experiments for >> Aqq^ at HERA 
and the nuclear shadowing experiments at Fermilab and CERN. For an excellent 
introduction to the field, see Ref. The results of the HERA experiments show 
a very rapid rise in parton distributions for x << 1 which is explained both by 
the conventional (operator product expansion), leading twist Double Leading Log 



approximation [y^,|l3| and the less conventional BFKL equation |15, i^- However, in 
the asymptotic limit of x — > 0, neither of these approximations are correct because 
they would both violate the unitarity bound on the growth of cross sections at 
asymptotic energies iQ. 

These explanations break down completely at very small x because the parton 
densities become very large and many body effects become important. Consequences 
of parton "overcrowding" are that two soft partons may recombine to form a harder 



parton or a parton may be screened by a cloud of surrounding wee partons |16, 17]. 
These processes inhibit the growth of parton distributions which saturate at some 
critical x. Indeed, these processes become important in nuclei at larger values of x 
than in nucleons. This may explain the strong A dependent shadowing seen in the 
deeply inelastic scattering (DIS) off nuclei at Fermilab and CERN fl^ . 

It may be argued that practically all one needs in order to determine the dynam- 
ics after a nuclear collision are the empirical nuclear structure functions at small 
X [jl9|. Using the QCD factorization theorem, products of these probabilities of 
finding a parton in the nucleus may then be convolved with the elementary parton- 
parton cross sections to determine parton scattering rates after the nuclear collision. 
However, factorization breaks down at small x (central rapidities) and coherence 
effects become important. Partons from one nucleus, at relevant transverse mo- 
mentum scales, do not resolve individual partons from the other nucleus. As in the 
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quantum theory of scattering, one needs to take the overlap of the wavefunctions- 
or more specifically, the small x Fock component of the nuclear wavefunction to 
determine the subsequent nuclear evolution. 

This question about the nuclear wavefunction is best formulated on the light 
cone using the method of light cone quantization One can write down the light 
cone QCD Hamiltonian, which is separable into a kinetic term and a potential term. 
Mueller has shown recently that for heavy quarkonia, where the scale of the coupling 
constant is set by the mass of the "onium" pair, light cone perturbation theory can 
be used to construct multi-parton eigenstates at small x [|^]. Note however, that 
despite many attempts, which go under the label "Light Front QCD", thus far the 
light cone approach to non-perturbative QCD has only had limited success |24|. 



We wish to argue that when the density of partons is extremely large, at very 
low X in a nucleon or in extremely large nuclei, the density of partons sets the scale 
for the running of the coupling constant. In other words, if 

then as{p) « 1. Here we will discuss specifically the application of weak coupling 
techniques in large nuclei at small values of x, with x « A^^^^. An intrinsic scale 
in the problem is set by the quantity fi^ ~ A^^^ fm~^, which is the valence quark 
color charge squared per unit area. Since it is the only scale in the problem, the 
coupling constant will run as a function of this scale |^ . 

In Section 2, we write down a partition function for the parton distributions at 
small X in the presence of the valence quarks which play the role of external sources 
in the problem. The background field for this theory is the non-Abelian analogue 
of the well known Weizsacker-Williams field in quantum electrodynamics |^5|. The 
parton distribution functions are formally expressed as correlation functions of a 2- 
dimensional Euclidean field theory with the effective coupling a^/i. The correlation 
functions to each order in as, involve an infinite resummation to all orders 'magfj, 0] . 



It was hoped that this classical theory would generate a screening mass (oc Os/u) 
which regulates the infrared behavior of the full theory. Recent lattice results 
suggest however that the classical problem may not be well defined in the infrared. A 
possible resolution is that the mass scale for the low momentum modes is generated 
at the quantum level by the high momentum modes ||2^. Higher order corrections 
to the background field are also discussed briefly in Section 2. 

Nuclear collisions are addressed in Section 3. Within the above picture, nuclear 
collisions can be understood as the collision of two Weizsacker- Williams fields. Since 
the fields are non-Abelian, the classical gluon field generated after the collision is 
obtained by solving the non-linear Yang-Mills equations with boundary conditions 
specified by the Weizsacker- Williams field of each nucleus In the central 

region of the collision, one therefore sees the highly non-perturbative (in agfi) evo- 
lution of the Weizsacker- Williams glue (and sea quarks). The time scale for the 
dissipation of these non-linearities is on the order of ~ I/q^/U. On time scales much 
larger than this time scale, the evolution of these fields can be described by the 
hydrodynamic scenario put forward by Bjorken. The quantum picture of nuclear 
collisions is discussed briefly, with particular reference to the "Onium" picture of 
Mueller. 

In Section 4, we will briefly discuss points of commonality as well as difference 
between the Weizsacker- Williams model and the above mentioned models as re- 
gards their conceptual foundations as well as their predictions for the experiments 
which will be performed at RHIC and LHC. These include the parton cascade model 
of Geiger and Muller p^,p9|], the various string fragmentation models ||3^ and the 
"color capacitor" models which describe particle production in ultrarelativis- 
tic nuclear collisions by the QCD analog of the Schwinger mechanism in quantum 
electrodynamics . 

Section 5 will contain our conclusions. 
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2 Computing parton distributions for a large nucleus 

In this section, the problem of calculating the distributions of partons in the nuclear 
wavefunction is formulated as a many body problem. We work in the infinite mo- 
mentum frame using the technique of light cone quantization. Our gauge of choice 
will be light cone gauge A'^ = 0. For an excellent discussion of the advantages of 



light cone quantization, we refer the reader to Ref. |2C|. In light cone quantization 
and light cone gauge, the electromagnetic form factor of the hadron F2, measured 
in deeply inelastic scattering experiments, is simply related to parton distributions 



by the relation [21 



F2{x, Q2) = ( p d%dk+xS{x - 1^) E 4«a) . (2) 



A=± 

We use light cone co-ordinates (x^ = (tibx)/\/2)- In the above, is the momen- 
tum of the nucleus, and kt are the parton longitudinal and transverse momenta 
respectively, x is the light cone momentum fraction, is the momentum transfer 
squared from the projectile and a^a is the number density of partons in momentum 
space. One only need integrate the calculated distributions up to the scale to 
make comparison with experiment. 

2.1 A partition function for wee partons in a large nucleus 

In QED, the infinite momentum frame wavefunction of the system with the external 
source in Eq. (||) is a coherent state Q]. Failing to do the same in QCD, we compute 
ground state expectation values instead. The partition function for the ground state 
of the low X partons in the presence of the valence quark external source is 

Z =< 0\e'^^~ \0 >= lim V < iV|e*^^" >q . (3) 

The sum above also includes a sum over the color labels of the sources of color 
charge (denoted by Q) generated by the valence quarks. 



The light cone Hamiltonian P (generator of translations in light cone time 
x^), in the presence of an external source, can be written as 

+ }-^\M- Pt)^{M+ Pt)^!;, (4) 

where P~ can be split into kinetic and potential pieces. Here p_p is the charge density 
due to the external source plus the dynamical quarks, At and -0 are dynamical 
components of the vector and spinor fields respectively, Et = d-At is the transverse 
electric field and Ff is the transverse field strength tensor. See Ref. Q for more 
details regarding the above expression and the conventions used. 

Valence quarks are predominantly found at large values of x. It is therefore 
reasonable to assume that they constitute the sources of the external charge seen 
by the wee partons. This current takes the form 

J^ = 5^'+Pa{x+,x^)5{x-). (5) 

In the gauge A'^ = 0, the static component is the only large component of 
the valence quark current. The transverse and minus components are proportional 
to and are therefore small. The current seen by the wee partons is pro- 

portional to (5(x~) if the valence quarks are Lorentz contracted to a size which 
is much smaller than a co-moving wee parton's wavelength. This is satisfied if 
2i?m/P+ << 1/xP^ =^ X « 1/Rm ~ A^^^^ where R is the nuclear radius and 
m the nucleon mass. 

Evaluating the trace in the partition function for quantized sources of color 
charge is difficult. We simplify the problem by resolving the transverse space as a 
grid of boxes of size cPxt » 1/pvai (or equivalently, parton transverse momenta 
Qf « Pvai) which contains a large number of valence quarks and hence a large 
number of color charges. This allows us to treat the sum over color configurations 
classically [Q. We average over the color charges by introducing in the path integral 
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representation of the partition function the Gaussian weight 

where p is the color charge density (per unit area) and the parameter fjp is the 
average color charge density squared (per unit area) in units of the coupling constant 
g. It can be written as 



2 

= Pval 



3A 4 



<Q">^'^2l9"- l.lA'/'fm-\ (7) 

where < >= 4(7^/3 is the average color charge squared of a quark. 

We can now write the partition function Z in the Light Cone gauge A_ = as 



Z = y [dAtdA+][d'tpUi;][dp] 

exp(^iS + ig j d'^xA^{x)5{x~)p{x) - J d^xtp'^{0,xt)^ . (8) 

Hence, the result of our manipulations is to introduce a dimensionful parameter 
p'^ ~ 1.1 A^^^ fm~^ in the theory. If we impose current conservation, D^J^ = and 
integrate over the external sources /), we obtain a non-local theory Q containing 
modified propagators and vertices. 

2.2 The classical background field of a nucleus 

Equivalently, one can find the classical background field in the presence of the 
external sources, compute correlation functions in this background field and finally 
integrate over the Gaussian random sources. 
The equations of motion are 

D^F^^'' = gr ■ = 5+''pa{x+,xt)5ix-) . (9) 

The background field which satisfies the above equation of motion is A'^ = 0, 
Ai{x) = 9{x~) ai{xt). Also, F12 = and V • q = gp{xt). Because the field strength 



Fi2 = 0, Ui is a pure gauge: r • = —j-UViU^. Combining the two equations 

''9 

results in the highly non-linear stochastic differential equation 

V • UVU^ = -ig^p{xt) . (10) 

Here p is the surface charge density associated with the current J . There is no 
dependence on x~ because we have factored out the delta function. The dependence 
on goes away because of the extended current conservation law D^J^^ = in the 
background field. 

To compute correlation functions associated with our classical solutions, we 
must solve the above equation and integrate the rho-dependent gauge fields over 
all color orientations of the external sheet of charge. In the matrix notation, 

(af (x,)af^'(0)) = j[dp]{u{x,)VU\x,)y^{u{Q)VU\G)y^^ 

X exp (^-^ j d^xtp''{xx)p''{xA_)^ . (11) 

The relation between distribution functions and propagators is straightforward 
and is discussed explicitly in Ref. |5|. 

where the trace is over both Lorentz and color indices. 

It was believed previously that the distribution function has the general form 

1 dN (iV,2-l) 1 1 ^.,2, 2 2^ l,o^ 

H{kjasfi), (13) 



7ri?^ dxd^kt vr^ x as 

where H{k1 / Ogp?) is a non-trivial function obtained by explicitly solving Eq. (11). 
It was further believed that in the "weak coupling" limit of kt » g'^p, H{y) 1/y 
and one obtains the Weizsacker-Williams result 

1 dN as/i2(iV2_i) 1 



(14) 
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scaled by fi"^. It was expected that the function H{k^ /a^/x^) would have in the strong 
coupling region the behavior a1fj.'^/{kf + M^) where Mg ~ a^/i is a screening mass 
which would regulate the divergence of the distribution function at small kt . If such 
a screening mass did exist, it would provide a simple understanding of saturation 
already at the classical level. 

Unfortunately, this turns out not to be the case. Unable to find an analytic 
solution to Eq. (10) we recently solved the correlation functions numerically on the 
lattice using the conjugate gradient method [^6|. Our preliminary results suggest 
the following. Weak (strong) coupling on the lattice holds when 0.2g^fiL << (>>)!, 
where L is the lattice size. In the weak coupling limit, our results indicate a discrete 
transverse momentum dependence which is of the 1/k^ Weizsacker-Williams form. 
As one increases g^fJ-L, there is an additional "transverse" contribution to the cor- 
relation function which displays an exponential fall off determined by a "screening 
mass". However, the amplitude of this term appears to diverge as ~ {g^^L)'^. At 
large values of g^^L, it appears that no solutions of the stochastic equations exist 
and the classical theory is ill defined in the infrared. 

Albeit it appears that no infrared stable screening mass is generated at the 
classical level, such a screening mass may still be generated by quantum fluctuations 
around the high transverse momentum modes Another possibility is that the 
Gaussian weight has a momentum dependence which modifies the infrared behavior 
of the correlation functions. Finally, the assumption that the source is a delta 
function on the light cone may be too severe. A way to deal with this problem is 
suggested in very interesting recent work by Balitsky [^3|. These ideas are outside 
the scope of the present work. We shall assume in the following that the computation 
of small fluctuations in the background field is not modified seriously by our lattice 
result (i.e., the theory can, for instance, be regularized by integrating over the 
random sources with an appropriate weight). 
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2.3 Quantum corrections to background field 



We now outline a procedure (the familiar Dyson-Schwinger expansion ||3^) to sys- 



tematically compute quantum corrections to our background field to all orders. The 
fully connected two point function is given by the relation 



In the above, {Ad) is the expectation value of the classical field to all orders in h. 

It can be expanded as {Ad) = A^^^ + where A^^^ is the solution discussed 

in Section 2.2. The term {AgAq) is the small fluctuation Green's function computed 
to each order in the classical field. The symbol < ■ ■ ■ >p indicates that we have to 
average over the external sources of color charge with the Gaussian weight described 
previously. We will briefly discuss the computation of the small fluctuation Green's 
function and how it may be used to compute the one loop correction to a) the gluon 
distribution function and b) the classical field A^^\ 

We begin by considering small fluctuations around our classical background, 
Ad = A^i + 5 A. Substituting this in the partition function in Eq. (|8|), we only keep 
terms 0(6^^) in the action. The small fluctuations propagator may be computed 
directly from the action or directly from the small fluctuation Yang-Mills equations, 
keeping terms linear in 5A, and solving the resulting eigenvalue equation. The final 
expression which involves several subtle features of light cone quantization is quite 
lengthy and the reader is referred to Ref. for the details. 

In Ref. we have used the small fluctuation Green's function in light cone 
gauge to compute the one loop correction to the background field as well as the one 
loop contribution to the classical field. The perturbative expression for the gluon 
distribution function to second order in is 




(15) 



ttW^ dxdPh 



1 dN 




1 + 



2asN^ 



vr 



In 




10 



Equation (16) contains both ln(l/x) and ln(A;t) corrections to the l/{xk^) dis- 
tribution and they represent the first order contributions to the perturbative ex- 
pansion for the distribution function. In the kinematical region of vahdity, these 
corrections are large. This signals that in order to properly account for the per- 
turbative corrections one has to devise a mechanism to isolate and sum up these 



leading contributions. This work is in progress |27|. 

The one loop corrections to the classical background field are computed as 
follows. We start with the classical equations of motion D^F^'^ = gj^ and expand 
the full gluon field as A'-'' = S^' + where is the background (classical) field, 
that is < > = B^^, while b^^ is the fluctuation (quantum) field with < 6'^ > = 0. 
Keeping up to quadratic terms in b^, one can write equations for the +, — and 
transverse components of the equations of motion. All the terms involving bilinear 
products of b^ in the minus and transverse components of the equations can be 
shown to vanish and these equations are identical to their classical counterparts. 
The equation for the -|- component is 

- d-d+B- - {Did+B')a = git + g<Jt>, (17) 



where the induced current in the above equation is related simply to the small 



fluctuations propagator GVix^y) by the equation 



d 

jt{x) = fabc < bl{x)d+bl{x) >= ifabc lim -^GUix, y) . (18) 

We find that the Fourier transform oi j^ix) is gj + (p)=n+:(p)Af)(p), where 
(p) is given by 

Ki{p) = 9Wr-^)Sab, (19) 



V 167r2 

which is the standard expression for the +i components of the polarization operator 



in light cone gauge |42|. We conclude from the above that the modifications to 



the background field introduced by the quantum fluctuations do not induce extra 
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terms in the expression for the distribution function [H]. This is consistent with the 



theorem of Dokshitzer, Diakonov and Troyan [38|. The effect of quantum corrections 
to the background field can be included by replacing the coupling constant g by the 
renormalized coupling constant gji which runs as a function of fj?. The structure of 
the background field at one loop remains unchanged. 

3 Nuclear collisions of Weizsacker— Williams fields 

In the previous section we discussed the properties of the Weizsacker- Williams field 
of a single nucleus. Recently, A. Kovner, L. McLerran and H. Weigert ||,10| have 



made significant progress in solving the classical problem of the evolution of these 
fields after the nuclear collision. 

Below, we outline very briefly their key results and refer the interested reader 
to their papers for further details. Before the two nuclei collide (for times t < 0), 
the Yang-Mills equations for the background field of two nuclei on the light cone is 
simply = and 

A' = 9{x-)e{-x+)a\{x±) + e{x+)9{-x-)ai{x±) (20) 

The two dimensional vector potentials are pure gauges (as in the single nucleus 
problem!) and for i < solve V • ai^2 = gpi,2ix±)- The interesting aspect of this 
solution is that the classical field configuration does not evolve in time for t < 0! 
This is a consequence of the highly coherent character of the wee parton clouds in 
the nuclei. 

The above solution for t < is a fairly straightforward deduction from the 
single nucleus case. What is very interesting is that the above mentioned authors 
find a solution to the field equations after the nuclear collision ( for t > 0). It is 
given by 

^± = ±x^a{T,x±) ; A' = a^{T,x±) , (21) 
12 



where r = \/t'^ — = \/2x+x^. The relation between fohows from the gauge 
condition x^A^ + x^A'^ = 0. This solution only depends on longitudinal boost 
invariant variable r and has no dependence on the space-time rapidity variable 
y = ^ In . This suggests that the parton distributions will be boost invariant for 
all later times. This result therefore justifies Bjorken's ansatz for the subsequent 
hydrodynamic evolution of the system. 

The above ansatz for the background field can be substituted in the Yang-Mills 
equations to obtain highly non-linear equations for a{T,x±) and a\{T,x±). The 



detailed expressions are given in Ref. [10|. The initial conditions for the evolution 



of these equations depend on the single nucleus solutions: 

a^\r=o = a{ + 4 ; a|r=o = y K,^] > (22) 

where a\ 2 are the background fields for the two nuclei. 

The Yang-Mills equations with the above boundary conditions are solved per- 
turbatively, order by order, by expanding the fields in powers of the valence quark 
charge density p. For asymptotically large r, Kovner et al. find that a gauge 
transform of the fields a and (denoted here by e and e\ respectively) have the 
form 



where oi and 02 can be expressed in terms of the p fields. In this equation, the 
frequency is a; =| k± \ and the vector = e^^k^ /to. The notation C. C. denotes 
complex conjugate. 

With the above form for the fields, the expressions for the parton number 



densities is straightforward. For late times, near z = 0, one obtains |41| 

dN 1 



3 



dy(Pk± (27r) 
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E I I' (24) 



Averaging over the p fields with the Gaussian weight in Eq. (6), one obtains the 
following result for the gluon distribution at late times after the nuclear collision: 

As suggested by the logarithm in the above equation, the transverse momentum 
integrals are infrared divergent. They are cut off by a mass scale a^/U which was be- 
lieved to result from the non-perturbative behavior of the classical fields. However, 
as we have discussed earlier, our recent lattice calculations suggest instead that the 
theory is not well defined in the infrared. Any such dynamically generated mass 
will therefore arise only at the quantum level. 

This brings up a related issue. Thus far we have only discussed the dynamical 
evolution of the classical fields. What about quantum effects? One way to include 
these is to do what we did for a single nucleus-look at small fluctuations around the 
background field of two nuclei The background field in this case is much more 
complicated than in the single nucleus case and the quantum problem is significantly 
more difficult. Another approach is to consider what quantum effects do to the 
coherence of the initial wavepacket. 



In this regard, A. H. Mueller's [22| formulation of the low x problem is relevant. 
He considers an "Onium" (heavy quark-anti-quark) state of mass M for which 
as{M) ^ 1. In weak coupling, the n- gluon component of the onium wavefunction 
obeys an integral equation whose kernel in the leading logarithmic and large Nc limit 



is precisely the BFKL kernel [15|. The derivation relies on a picture in which the 
onium state produces a cascade of soft gluons strongly ordered in their longitudinal 
momentum; the i-th emitted gluon has a longitudinal momentum much smaller 
than the i — 1-th. 

In the large limit, the n gluons can be represented as a collection of n-dipoles. 
Hence, in high energy onium-onium scattering, the cross section is proportional to 
the product of the number of dipoles in each onium state times the dipole-dipole 
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scattering cross section [^]. This cross section is given by two gluon exchange (the 
pomeron). More comphcated exchanges involving multi-pomeron exchange have 



been studied recently by Salam |3C]. However, despite the mathematical elegance 
and simple interpretation of the onium approach, it is unclear whether it can be 
extended to nuclei. 

4 Parton cascades and color capacitors 

In this section we will briefly discuss, in relation to the model discussed in earlier 
sections, some other attempts to model the initial conditions for ultrarelativistic 
heavy ion collisions. They may be broadly (and somewhat imprecisely) classified as 
follows: a) perturbative QCD based models which assume the factorization theorem 
and incoherent multiple scattering to construct a spacetime picture of the nuclear 
collision, and b) non-perturbative models where particle production is based on 
string fragmentation or pair creation in strong color fields. 

Among perturbative QCD based models, the parton cascade model of Geiger 



and Miiller |28,2£] has been applied extensively to study various features of heavy 
ion collisions. The evolution of classical phase space distributions of the partons is 
specified by a transport equation of the form 



dt dr 



Fa{p,r,t) = Ca{p,f,t), (26) 



where Fa are the classical phase space distributions for particle type a and Ca is 
the corresponding collision integral. The matrix elements in the collision integral 
are computed from the relevant tree level diagrams in perturbative QCD. 

The initial conditions in the parton cascade model are specified at some initial 
time t = to by the distribution Fa{p, f,t = to) = Pa{p, P)Ra{r, R)- The momentum 
distribution Pa{p, P) = fa(x,QQ)gipt) is decomposed into an uncorrelated product 
of longitudinal and transverse momentum distributions respectively, where fa{x, Qq) 
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is the nuclear parton distribution which is taken from Deep Inelastic Scattering 
experiments on nuclei at the relevant Qq and g{pt) is parametrized by a Gaussian 
fit to proton-proton scattering data. The spatial distribution of the partons is 
described by a convolution of a Woods-Saxon distribution of nucleons in the nucleus 
and an exponential distribution of individual partons within each nucleon. Details 
regarding both initial conditions may be found in Ref. ||2^. Another model which 
takes as input the perturbative QCD cross sections is the HIJING model pS], pO|] 
which describes nuclear scattering in an eikonal formalism which convolves binary 
nucleon collisions. In both models, detailed predictions have been made for various 
observables at RHIC-in particular, for mini-jet production. 

As suggested by the above, both models make assumptions which are not nec- 
essarily motivated by perturbative QCD. To an extent, this is inevitable because 
one is forced to model the soft physics. Where these approaches differ significantly 
from the Weizsacker-Williams approach is in the factorization assumption, namely, 
that partons from one nucleus resolve individual partons of the other in each hard 
scattering. We have argued that the small x partons which dominate the physics 
of the central region instead have highly coherent "wave- like" interactions. This 
results in a vastly different space-time picture for the nuclear collision-at least for 
the very primordial stage of the nuclear collision. 

Naturally, the predictions of these models will differ significantly from the 
Weizsacker-Williams model. For instance, because of the intrinsic pt ^ ^ carried 
by the Weizsacker-Williams (or "equivalent") gluons, gluon production is enhanced 
by a factor ag relative to the lowest order gg ^ gg mini-jet process in a cascade. A 
simple explanation for this enhancement is that because the valence quarks absorb 
the "recoil", two off-shell equivalent gluons can combine to produce an on-shell 
gluon. This will impact significantly the many signatures to be studied at RHIC 
and LHC such as jet production and dilepton and photon production. Further, the 
intrinsic pt of the gluons ensures that intrinsic charm and strangeness production is 
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significantly larger in the Weizsacker- Williams model 

The non-perturbative models [^0| primarily attempt to describe the soft physics 
in ultrarelativistic nuclear collisions so it is not clear that there is much overlap with 
the Weizsacker-Williams model. However, the latter does provide some insight into 
one of these approaches, which we shall dub the "color capacitor" approach. Here 
it is assumed that the nuclei produce a homogeneous chromo-electric field which 
produces particles non-perturbatively by a mechanism analogous to the Schwinger 
mechanism for strong electromagnetic fields. The evolution of these fields (including 
back-reaction) is determined by a Boltzmann-like equation where the source term 
now is given by the pair production rate pl|,p^. 

An important assumption in these color capacitor models is that of homogene- 
ity of the initial field configurations. However, the results discussed in the previous 
section suggest that the Yang-Mills fields are highly non-linear and inhomoge- 
neous. The time scale r >> l/a^/i is the time scale for the dissipation of the 
non-linearities in the fields. It would be interesting to see how the solutions to the 
transport equations are modified for initial conditions given by the inhomogeneous 
Weizsacker-Williams field configurations. 

5 Conclusions 

We have described in this Comment a QCD based approach to describe the ini- 
tial conditions for ultrarelativistic nuclear collisions. The central region of these 
collisions is dominated by "wee" partons which carry only a small fraction of the 
nuclear momentum. We have argued that for very large nuclei these partons are 
only weakly coupled to each other. However, due to their large density, many body 
effects are extremely important. The classical behavior of these quanta (the QCD 
analogue of the Weizsacker-Williams equivalent photons) can be described by an 
effective two dimensional field theory. Quantum effects are treated by constructing 
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the small fluctuations propagator in the background field of these quanta. 

An important objective of this approach is to understand if there is a "Lipatov 
region" in nuclei where the parton densities grow rapidly and if the shadowing of 
parton distributions in nuclei can be understood to result from the precocious onset 
of parton screening. It is probable that deep inelastic scattering experiments off 
large nuclei will be performed at HERA in the near future [Q. If so, one may ex- 
pect unprecedentedly high parton densities and interesting and perhaps unexpected 
phenomena in these experiments. 

These DIS experiments on nuclei at HERA would nicely complement the heavy 
ion program at RHIC and especially LHC since they probe the same range of Bjorken 
X. The results of these experiments would therefore place strong bounds on mini- 
jet multiplicities and other signatures of nuclear collisions. Note that these observ- 
ables are extremely sensitive to the initial parton distributions (for a discussion, 
see Ref. [^]). However, to fully understand the dynamics of nuclear collisions at 
central rapidities, we have to understand the initial conditions ab initio-preferably 
in a QCD based approach like the one discussed in this paper. 

At the moment there are still many open questions which remain unresolved. 
An empirical question is with regard to the applicability of weak coupling methods 
to large nuclei. Is the bare parameter fi^ ~ A^^^ fm"^ large enough? One may 
argue on the basis of Renormalization Group arguments that this parameter should 
effectively be larger and should grow with the increasing parton density at small x. 
However, these arguments are not rigorous at this stage. 

A more serious problem is suggested by recent lattice simulations of the 2- 
D effective field theory which show that the classical correlation functions diverge 
quadratically with the lattice size L. Identifying why this divergence occurs and 
how the background field may be modified accordingly needs to be resolved satis- 
factorily. Finally, the problem of Lipatov enhancement and saturation in nuclear 
parton distributions is not yet settled. 
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Despite the many technical problems that remain, there is much cause for opti- 
mism since it appears now that the problem of initial conditions in ultrarelativistic 
nuclear collisions can be treated systematically in a QCD based approach. Because 
the various empirical signatures depend sensitively on the initial conditions, one 
may hope to identify and interpret the elusive quark gluon plasma in ultrarelativis- 
tic nuclear collisions at RHIC and LHC early in the next millenium. 

Acknowledgments 

Foremost, I would like to acknowledge my debt to my collaborators A. Ayala, R. 
Gavai, J. Jalilian-Marian, A. Kovncr, H. Weigert and especially L. McLcrran. I 
would also like to acknowledge useful discussions with I. Balitsky, B. Miillcr, R. 
Narayanan, J. Qiu, T. Schaefer and M. Strikman. Finally, I would like to thank J. 
Kapusta for his patience and good humor. 

Research supported by the U.S. Department of Energy under grants No. No. 
DOE Nuclear DE-FG06-90ER-40561. 

References 

[1] H. Cheng and T. T. Wu, Phys. Rev. Lett. 24 759 (1970). 

[2] see for instance the proceedings of Quark Matter 95, Ed. A. M. Poskanzer, J. 
W. Harris and L. S. Schroeder, North Holland Publishers. 

[3] L. McLerran and R. Venugopalan, Phys. Rev. D49 2233 (1994). 

[4] L. McLerran and R. Venugopalan, Phys. Rev. D49 3352 (1994). 

[5] L. McLerran and R. Venugopalan, Phys. Rev. D50 2225 (1994). 

[6] A. Ayala, J. Jalilian-Marian, L. McLerran and R. Venugopalan, Phys. Rev. 
D52 2935 (1995) . 

19 



[7] A. Ayala, J. Jalilian-Marian, L. McLerran and R. Venugopalan, Phys. Rev. 
D53 458 (1996) . 

[8] R. Venugopalan, Nucl. Phys. A590 147c (1995). 

[9] A. Kovner, L McLerran and H. Weigert, Phys. Rev. D52 6231 (1995). 

[10] A. Kovner, L McLerran and H. Weigert, Phys. Rev. D52 3809 (1995). 

[11] B. Badelek, K. Charchula, M. Krawczyk and J. Kwiecinski, Rev. of Mod. Phys. 
64 927 (1992). 

[12] D. J. Gross, Phys. Rev. Lett. 32 1071 (1974). 

[13] R. D. Ball and S. Forte, Phys. Lett. B351 313 (1995). 

[14] M. Proissart, Phys. Rev. 123 1053 (1961). 

[15] E.A. Kuraev, L.N. Lipatov and Y.S. Fadin, Zh. Eksp. Teor. Fiz 72, 3 (1977) 
{Sov. Phys. JETP45, 1 (1977) ); LA. Balitsky and L.N. Lipatov, Sov. J. Nucl. 
Phys. 28 (1978) 822. 

[16] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 1 (1983); E. Laenen 
and E. M. Levin, Nucl. Phys. B451 207 (1995). 

[17] A. H. Mueller and J. Qiu, Nucl. Phys. B268 427 (1986). 

[18] T. Carroll, A. Witzmann and H. Melanson, in Proc. XXVIII Rencontres de 
Moriond, Les Arcs, Savoie, France 1993. 

[19] K. Eskola, Nucl. Phys. A590 383c (1995). 

[20] S. J. Brodsky and H. -C. Pauli, Lectures at 30th Schladming Winter School in 
Particle Physics, Schladming, Austria, March 1991. 



20 



[21] Yu. L. Dokshitzer, V. A. Khoze, A. H. Mueller and S. I. Troyan, Basics of 
Perturbative QCD, Editions Frontiers, 1991. 



[22; 
[23 
[24 
[25; 

[26 
[27; 

[28 
[29 
[30 



[31 

[32; 

[33 
[34 



A. H. Mueller, Nucl. Phys. B415 373 (1994); Nucl. Phys. 437 107 (1995). 
A. H. Mueller and B. Patel, Nucl Phys. B425 471 (1994). 
M. Burkardt, submitted to Adv. Nucl. Phys., [tiep~ph/9505259 



J. D. Jackson, Classical Electrodynamics, Section 15.4, John Wiley Publishing 
Co. 

R. Gavai and R. Venugopalan, to be published. 

A. Kovner, L. McLerran, H. Weigert and J. Jalilian-Marian, private commu- 
nication. 

K. Geiger, B. Miiller, Nucl. Phys. B369 600 (1992). 
K. Geiger, Phys. Rep. 258 237 (1995). 

A. Capella, U. Sukhatme, C-I. Tan and J. Tran Thanh Van, Phys. Rep. 236 
(1994); B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Phys. Rep. 
97 31 (1983). 

J. M. Eisenberg, LBL-34831, Proceedings of Workshop on Pre-Equilihrium 
Parton Dynamics, Berkeley, August 23-September 3, 1993. 

F. Cooper, J. M. Eisenberg, Y. Kluger, E. Mottola, B. Svetitsky, Phys. Rev. D 
48 190 (1993). 

I. Balitsky, preprint MIT-CTP~2470, September 1995. 
J. D. Bjorken, Phys. Rev. D27 140 (1983). 



21 



[35] L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics, Addison- Wesley 
Publishing Co., Inc. 1989. 



[36 
[37; 
[38 
[39 

[4o; 

[41 
[42; 
[43' 
[44; 



G. P. Salam, preprint Cavendisli-HEP-95/05, |hep-ph/9509353 
A. Makhlin, Phys. Rev. C52 995 (1995). 

Yu. Dokshitzer, D. Diakonov and S. Troyan, Phys. Rep. 58 270 (1980). 

K. Eskola and Xin-nian Wang, Phys. Rev. D49 1284 (1994). 

M. Gyulassy and Xin-nian Wang, Comput. Phys. Commun. 83 307 (1994). 

J.-P. Blaizot and A. H. Mueller, Nucl. Phys. B289 847 (1987). 

G. Leibbrandt, Rev. of Mod. Phys. 59 1067 (1987). 

V. Del Duca, DESY preprint, |hep-ph/9503"22^ . 



For a discussion, see proceedings of the workshop on Nuclei at HERA and Heavy 
Ion physics, Ed. S. Gavin and M. Strikman, Brookhaven, november 1995. 



22 



